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ABSTRACT 

We present broadband (3 - 78 keV) NuSTAR X-ray imaging and spectroscopy of the Crab nebula 
and pulsar. We show that while the phase-averaged and spatially integrated nebula + pulsar spectrum 
is a power-law in this energy band, spatially resolved spectroscopy of the nebula finds a break at ^9 
keV in the spectral photon index of the torus structure with a steepening characterized by AT ^ 0.25. 
We also confirm a previously reported steepening in the pulsed spectrum, and quantify it with a 
broken power-law with break energy at ^12 keV and AT ~ 0.27. We present spectral maps of the 
inner 100 " of the remnant and measure the size of the nebula as a function of energy in seven bands. 
These results find that the rate of shrinkage with energy of the torus size can be fitted by a power-law 
with an index of 7 = 0.094 ± 0.018, consistent with the predictions of Kennel and Coroniti (1984). 
The change in size is more rapid in the NW direction, coinciding with the counter-jet where we find 
the index to be a factor of two larger. NuSTAR observed the Crab during the latter part of a 7 -ray 
flare, but found no increase in fiux in the 3-78 keV energy band. 

Subject headings: SNR: PWN — Crab - X-rays: ISM - stars: neutron - pulsars: individual (Crab) 


1 . INTRODUCTION 

The Crab is the prototypical Pulsar Wind Nebula 
(PWN), characterized by a center-filled synchrotron neb¬ 
ula that is powered by a magnetized wind of charged par¬ 
ticles emanating from a central ly located pulsar formed 
during the supernova explosion (jWeiler & Panagia|1978 ). 
Due to its brightness, proximity of ^2 kpc, and well- 
known explosion date of 1054 AD, the Crab is the best 
studied object of its kind, and the literature is rich with 
hundreds of publications and reviews describing its p rop¬ 
erties across all energy ban ds (see e.g. Hester (2008) and 
Biihler & Blandford (2014) for a recent review)"! 


Detailed ima gesTrom the Hubble Space Telescope (Hes- 
ter et al.|1995 ) and Chandra ( Weisskopf et al.| 2000 ) have 
revealed the nebula’s morphological complexities. In the 
optical band the remnant measures ^3' across its longest 
axis with thermal filaments composed of ejecta from the 
explosion confining the synchrotron nebula. In X-rays 
the remnant is considerably smaller and shows both torus 
and jet structures. The symmetry axis is tilted at about 
27° to the plane of the sky with the NW edge closer to the 
observer (Ng & Romani 2004). The jet emerges toward 
the observer to the bE, and a less collim ated structure 
(the co unter-jet) extends away to the NW. Ng & Romani 
(2004) fit simple models to the morphology to obtain a 
radial flow velocity through the torus of 0.550 ± 0.001c, 
and at these velocities relativistic beaming brightens the 
nearside (NW) of the torus. 


The wind likely te rminates at a shock z one about 
10 " from the pulsar (| Weisskopf et al.|| 2000 |), and the 
post shock particles and magnetic fields propagate non- 
relati vistically eithe r by diffusion (|Gratton||1972| iWilson 
1972| or ady ection (|Rees Cunn|l974[ |Kennel & Coro- 
niti||l984a|b[) to the edge of the remnant, emitting their 


energy as synchrotron radiation. Both the process of par¬ 
ticle acceleration in relativistic shocks and the transport 
of particles and fields downstream, occur in other astro- 
physical settings such as gamma-ray bursts and jets in 
active galaxies, and studying them in a relatively nearby 
spatially resolved source, may have application beyond 
the understanding of the Crab and PWNe in general. 

From radio to TeV, the emission from the Crab neb¬ 
ula is non-thermal and peaks in the optical through X- 
ray. At radio wavelengths the nebula’s integrated emis¬ 
sion is a power-law s pectrum with ind ex {Sj^ oc 
a = —0.299 ± 0.009 (Baars et al. 1977). In the opti¬ 


cal, the synchotron spectrum is steeper with a gradual 
turnover occurring so mewhere between 10 and 1000 /am 
( Marsden et al.|p^M ). In the X-ray band the Crab neb- 
ula-hpulsar phc)ton index is F ^ 2.1 and further soft - 
ens above 100 keV to ^2.23 (jJourdain & Roques||2009|). 
The Crab nebula component alone l ikewise softens above 
100keV to ^2.14 up to 300keV (Pravdo, Angelini & 


Harding 1997) and soft ens further to 2.22 7 ± 0.013 be¬ 
tween 0.75 and 30 MeV ( Kniper et al.| 2001 ) where it con¬ 


tinues to soften until 7001VleV, beyond which the inverse 
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TABLE 1 
Observations Log 


Obs ID 

Date 

l^xposure 

seconds 

off-axis angle 
arcminutes 

Comment 

1001302200^ 

2012-09-20 

2592 

1.5/2.0 

Crab22 

10013022004^ 

2012-09-21 

2347 

1.5/2.0 

Crab22 

10013022006^ 

2012-09-21 

2587 

1.5/2.0 

Crab22 

10013031002^ 

2012-10-25 

2507 

2.0/2.6 

Crab31 

10013032002^ 

2012-11-04 

2595 

1.2/1.3 

Crab32 

10013034002 

2013-02-14 

988 

0.8/1.5 

Crab34 

10013034004 

2013-02-14 

5720 

0.7/0.9 

Crab34 

10013034005 

2013-02-15 

5968 

0.9/0.6 

Crab34 

10013037002 

2013-04-03 

2679 

2.0/2.4 

Crab37 

10013037004 

2013-04-04 

2796 

2.9/3.5 

Crab37 

10013037006 

2013-04-05 

2944 

2.9/3.5 

Crab37 

10013037008 

2013-04-18 

2814 

3.0/3.7 

Crab37 

80001022002^ 

2013-03-09 

3917 

1.5/1.8 

CrabToo 

10002001002^ 

2013-09-02 

2608 

1.72/2.09 

CrabSci 

10002001004^ 

2013-09-03 

2386 

1.73/2.26 

CrabSci 

10002001006^ 

2013-11-11 

14260 

1.08/1.28 

CrabSci 


^Effective exposure time corrected for dead-time. 
^Subset used for deconvolution analysis. 


Comp ton component sets in (Meyer, Horns & Zechlin 


2010). The pulsed spectrum follows a ditterent spectral 


energy distrib ution (SEP) with a more complex spec¬ 
tral evolution ( l^iper et aL]|200 ^ [Weisskopf et al.||2Qll 


Kirsch et ah ||26o6| r' 


in Ge V 7 -rays the Crab has been observed to flare 
rapidly (Tavani et ah 2011 lAbdo et ah 2011) roughly 
once a year tor a duration or ^10 days, with a flux in¬ 
crease of a factor of ^5. The origins of these flaring 
episodes are currently not understood, but due to the 
rapidity of the flares and the fast cooling time of syn¬ 
chrotron radiation, the flare emission is presumed to be 
of synchrotron origin rather than inverse Compton or 
bremsstrahlung (Abdo et al.||2011 ). 

No single model has yet successfully been devised to 
explain the properties of the nebula across all energy 
bands. The ratio of wind flow times to the particle ra¬ 
diative lifetimes is a strong model dependent parameter 
tied to the mechanism of energy transport and is observ¬ 
able as an energy-dependent nebular size. The predic¬ 
tion of both diffusion and advection models is a decrease 
in size of the nebula with increasing photon energy due 
to higher-energy electrons dying out sooner than lower- 
energy ones. In X-rays the expectation is for a diffusion- 
driven nebul a to be smaller t han one that is advection 
dominated. Ku et al. (1976) first collected broadband 
evidence confirming shrinkage of the Crab nebula, and 
this limited data set pointed towards a combination of 
diffusive and advective electron transport. 

The spatial dependence of the broad-band spectrum 
provides another observation al constraint on m odels. De¬ 
tailed Chandra observations (Mori et al.|2004) show that 
the torus spectrum is roughly uniform with a steepening 


at the edge, in lin e with predictions made by [Kennel fc 
Coroniti ( 1984a|b ), and indicates that for X-ray-emitting 


particles, the transport in the Crab nebula appears to be 
dominated by advection rather than diffusion. 

Above 10 keV, observations with the ability to spa¬ 
tially resolve the Crab nebula have been limited to one¬ 
dimensional scanning techniques and lunar occultation, 
which have limited signal to noise and lack true two¬ 


dimensional imaging (jPelling et al.||1987|). The NuSTAR 
high-energy X-ray focusing mission has the ability to 
make sensitive imaging observations of diffuse sources 
above 10 keV for the first time. NuSTAR^s two co- 
aligned telescopes operate from 3-78 keV, and have an 
imaging r esolution of 18^^ FWH M and ^ V Half Power 
Diameter (Harrison et al.|2013). NuSTAR can image the 
Crab above 10 keV with sufticient resolution to investi¬ 
gate the spe ctral and spatial properties of the nebula. 

In Section [3.1 j we present the g lobal properties of the 
pulsar and nebula. In Section [3^ we discuss the NuSTAR 
observati ons d uring the y-ray haring on the 9th of March 
2013. In |3.3| we present phase resolved spectroscopy of 
the Crab in the energy range 3-78 keV, tracing the 
spectru m o f the pulsar as a function of its phase. In 
Section 3.4 we present spectral maps of the inner 100" of 
the nebula, using pulse-off phase intervals, and in Section 
|3.5| we present deconvolved NuSTAR images to investi¬ 
gate the physical size of the Crab as a function of energy. 
In the discussion section we review and summarize the 
new findings. 

2. OBSERVATIONS AND DATA REDUCTION 

NuSTAR has two co-aligned telescopes with corre¬ 
sponding focal planes EPMA and EPMB. Each focal 
plane consists of four hybrid hard X-ray pixel detectors 
with a total field of view of 12.5'x 12.5b The detectors 
are denoted by numbers 0 through 3, with the optical 
axis placed in the inner corner of detector 0. NuSTAR 
observed the Crab multiple times throughout 2012 and 
2013 as part of the instrument calibration campaign. The 
observations span a wide range of off-axis angles, some of 
which were deemed too far off-axis for the analysis pre¬ 
sented here, and Table lists the chosen subset. All of 
these observations have the pulsar within 4' of the opti¬ 
cal axis and located on detector 0, with a total elapsed 
live-time corrected exposure of 59.7 ks. The accuracy of 
the absolute source location can vary due to the thermal 
environment with 3a offsets typically on the order of 8" 
(Harrison et al. 2013), but because of the very peaked 
pulsar emission it is straightforward to register all obser¬ 
vations to the pulsar location. 

Despite being a bright target ( r\j 250 cts s ^ per EPM 
with a dead time of ^ 50%), no special processing or pile- 
up corrections are needed since the focal plane detectors 
have a triggered readout with a short (1 /is) preampli¬ 
fier shaping time Harrison et al. (2013|. The data were 
reduced using the NubTARDAb pipeline version vl.3.0 
and CALDB version 20131223 with all standard settings. 
We extracted spectra using standard scripts provided by 
the NuSTARDAS pipeline. 

The background for the Crab only becomes impor¬ 
tant at energies above 60 keV where the internal detector 
background dominates over sky background. Since the 
internal background varies from detector to detector, ide¬ 
ally one would extract background from the same detec¬ 
tor as the source. Unfortunately the brightness and ex¬ 
tent of the Crab precludes extracting a local background 
from the same detector. We therefore simulated the 


backgrounds using the nuskybgd tool (Wik et al. 2014) 
and created a master background for each observation 
for a 200" extraction radius. To test for fluctuations in 
the background, we ran numerous realizations fitting a 
typical representative spectrum. We found that the fits 
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Fig. 1. — Ratio of the data to the best-fit model for all the Crab 
observations listed in Table The structured residuals, such as 
those around 20 - 25keV are related to the calibration process 
(piece-wise linear spline interpolation), but are small (~2%), typi¬ 
cally broad and at energies not to be mistaken for emission lines. 
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Fig. 2.— Spectra and fit of FPMA assuming a power-law index 
of F = 2.1. The normalizations of the curves have been scaled 
to illustrates how the index of the Crab softens with increasing 
extraction radius until at an extraction radius of 200" the power- 
law index of 2.1 is recovered. 


were not sensitive to the background fluctuations. 

In the text and tables all errors are reported at 90% 
confidence, and all fits are perfo rmed with XSPECQ us¬ 
ing Cash statistics (Cash||1979) unless otherwise stated. 
Spectra shown are re binned tor display purposes only. 

3. DATA ANALYSIS 
3.1. Phase Averaged Speetroseopy 

The spatially integrated spectrum of the Crab nebula 
+ pulsar in the 1 - 100 keV X-ray band has been well- 
described by a power-law with photon index T ^2.1 
(RXTE, BeppoSAX , EXOSAT, INTEGRAL/:SEMAX) 
^Kirsch et al.(2005|). Above 100keV the hard X-ray in- 
struments [IJy TEURAL/SFl/lSGRl, CGRO) measure a 
softer index of T ^2.20 - 2.25, and below lOkeV in¬ 
struments with CCD detectors a harder spectrum. This 


^ http://heasarc.gsfc.nasa.gov/xanadu/xspec/ 


hardening in CCD X-ray instruments comes about from 
photon pile-up, and although models exist for these in¬ 
struments to deal with the pile-up, the Crab usually still 
challenges the models and requires special non-standard 
reductions. In addition, it is common practice for these 
instruments to excise the piled-up regions, which removes 
part of the integrated spectrum and thus can techni¬ 
cally no longer be directly compared to other instruments 
where this excision has not occurred. The non piled-up 
instruments covering the 1 - 100 keV band agree on a 
photon index T = 2.1 ± 0.02 and have not measured any 
curvature in the Crab spectrum across this band. 

Over the 16 years RXTE was operational and regularly 
monitoring the Crab, the spectral index was se en to vary 
by a peak-to-pea k variation of AT ^ 0.025 (Shaposh- 


nikov et al. |2012), perhaps due to magn e tosom c waves 
in the nebula (e.g., Spitkovsky & Arons (2004|)). This 
variation is consistent with the observed spread between 
instruments, but is slow and on average over the 16 years 
the Crab remained at T = 2.1. Because the average in¬ 
dex covers several instruments and the deviations from 
it are small, NuSTAR has calibrated the effective area 
against a Crab index of T = 2.1. 

A total of 39 Crab observations, spanning off-axis an¬ 
gles from 0-7', went into adjusting the effective area an¬ 
cillary response files (ART) of NuSTAR. This was done 
using a piece-wise linear spline interpolation as a func¬ 
tion of energy and off-axis angle. Cross-calibration cam¬ 
paigns on 3C273 and PKS2155-304 have been used to 
confirm that the ART adjustments have not introduced 
a systematic offset, and well-known power-law sources, 
such as CenA, to confirm they still appear as power-law 
sources. We used the quasar 3C 273 to calibrate the Nh 
column and derived a value of Nh = (2±2) xlO^^ cm“^ 
for the Crab. At 3keV NuSTAR is only very marginally 
sensitive to Nh columns of 10^^ cm“^, which is the rea¬ 
son behind the large error on the Nh column. We have 
frozen the column to the above for all fits, emplo ying 
Wilms abundances d Wilms, Allen fc McC ray||2QQQ|) and 
Verner cross-sections (jVerner et al.||l996j). An extensive 
discussion on the choice of N h and efiective area calibra¬ 


tion can be found in Madsen 


(prep). 


The Crab flux has been registered in multipl e instru 


ments to decline ov er a 2 year period of ^ 7% (Wilson- 


Hodge et al. 2011) across the 15 - 50keV banal Tins 


corresponds to a decay in the flux of 3.5% per year over 
the period it has been observed. We set the Crab nor¬ 
malization to 8.5 to optimize (and minimize) the cross¬ 
calibration constants between concurrent X-ray observa¬ 
tories {Ghandra, Swift, Suzaku, XMM-Newton). In the 
followed, however, the absolute flux or a variation of it, 
has no influence on the results. 

We combine all the observations from Table ^ exclud¬ 
ing those denoted by “Crab34” that have the central 
region of the source falling on the gap between detec¬ 
tors in a way that complicates flux correction. Figure [T 
shows the ratio of data to the best-fit power-law modd 
r = 2.0963 ± 0.0004 and cross normalization between 
FPMA and FPMB of 1.002. The structured residuals, 
such as those around 20 - 25 keV are related to the cal¬ 
ibration process and are small (^2%), typically broad 
and at energies not to be mistaken for emission lines. 

It is important to note that in the subsequent analysis, 
we are focusing on changes relative to the above mea- 
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sured spectrum. The data quality is high, and we can 
measure slight changes within the spectrum with great 
accuracy. This is not to be mistaken for our knowledge 
of the absolute value of the spectral index. Individual 
fits to the 39 Crab observations yields an average mea¬ 
sured photon-index of T = 2.1 with a la spread of ±0.01. 
However, this error on our absolute measurement is not 
relevant to our analysis since we are only concerned with 
relative changes. 

To illustrate the range of spectral change in the Crab’s 
integrated spectrum as a function of radius, we extract 
all events within circles of radius: 12", 25", 50", 75", 
150" and 200" centered on the Crab pulsar and subtract 
the simulated master background. Figure shows the 
ratio of these spectra to the canonical power-law model 
T = 2.1, where the relative normalizations have been 
set equal at 3 keV. As already discussed, this progressive 
softening of the index is predicted by theory and has 
been measured at all energies (up to lOkeV) where it is 
possible to spatially resolve the Crab. 


3.2. Crab flaring 

In addition to the long time scale flux variations ob¬ 
served in the X-ray band, the Crab i s known to flare in 
GeV 7 -rays. Observations wit h Agile (Tavani et al.|2011 ) 
and Fermi (|Abdo et al.||2011|) in the U.l - 1 GeV range 
have observed hares with a frequency of roughly ^ 1 a 
year with typical durations of ^ 10 days and a flux in¬ 
crease of about a factor of 5. The radiation is thought 
to be of synchrotron origin due to the rapidity of the 
flares and the fast electron cooling time, as opposed to 
Bremsstrahlung or inverse-Compton emission tha t have 
cooling tim escales of the order of ^ 10^ —10^ years (Abdo 


et al. 2011). The emission region is furt her thought to 
be Do ppler boosted towards the observer (Buehler et al. 
2012) and causality arguments suggest the flares orig- 


mate from a very small region, most likely inside the 
termination shock zone. However, the spatial resolution 
of current y-ray instruments is not sufficient to r esolve 
the inner parts of the nebula. |Weisskopf et al.| (|2013|) 
analyzed Chandra^ Keck and VLA data during the 2011 
April flare to look for such a counterpart, but none of 
these instruments found conclusive evidence of a change 
in the Crab emission. 

On the 9th of March 2013 NuSTAR tri ggered an ob¬ 
serva tion during one such flaring episode (Mayer et al. 
2Q13[ ). This observation is labeled ‘CrabToO' in Table 
T] and had a duration of ^16 ksec and an effective ex¬ 
posure of ^4 ksec after taking dead time, SAA passages, 
and occultation into account. NuSTAR caught the Crab 
at the tail end of flaring and did not detect any spectral 
variation or flux change beyond the level that can be ex¬ 
pected from calibration (±5% in flux and 0.01 in spectral 
index change). 


3.3. Phase Resolved Speetroseopy 
3.3.1. Pulsar Spectrum 

The Crab pulsar has a 33 ms period and a pulse profile 
that is double peaked in all energy bands. In the X- 
rays the primary peak is higher than the secondary, and 
the off-pulse period spans about 30% of the phase (see 
Figure 3, middle panel). The Crab ephemeris is routinely 



Fig. 3.— Top: The pulse profile (P = 33ms) of module A is shown 
in solid, the live-time curve is shown dashed, and the boundaries 
of the extracted phase bins as vertical dotted lines. The minimum 
of the live-time curve occurs just after the peak and the oscillatory 
pattern is due to the dead-time of 2.5 ms. Middle: Live-time cor¬ 
rected pulse profile. Bottom: Power-law index averaged between 
17 " and 50'' regions are shown in diamonds for Pi and crosses r 2 . 
The histogram shows the relative normalization of pulsed emission 
to the nebular emission (A^pulse/^nebula) within each bin. 


calculated by the Jodrell Bank Observatory and we 
used the closest ephemeris entry for each observation to 
obtain the pulse profile. 

For the phase resolved spectroscopy, all observations 
listed in Table were combined for FPMA, and for 
FPMB “Crab34” was excluded due to proximity of the 
pulsar to the chip gap. The Ancillary Response Files 
(ARFs) of the different observations were combined, 
while the same detector response matrix (RMF) could 
be used since the pulsar was on detector 0 for all obser¬ 
vations. We barycenter corrected event times using the 
barycorr routine that is part of the HEASOFT FTOOL 
library, and corrected the NuSTAR clock for thermal 
drifts. 

The resulting raw (not live-time corrected) pulse pro¬ 
file is shown in the top panel of Figure]^ (solid line). The 
pulse profile exhibits a distinct oscillatory pattern due to 

^ (http://www.jb.man.ac.uk/pulsar/crab.html) 
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dead time effects associate d with the 2.5 msec event read¬ 
out (jHarrison et al.||2013|). The decreased probability of 
detecting an event just after pulse peak results in a pat¬ 
tern of damped oscillations with a period of ^3.4 msec 
and causes the live-time fraction to vary throughout the 
pulse period as shown by the dashed curve. To cor¬ 
rect for this effect we use the “PRIOR” column in the 
event list, which specifies the elapsed time since the prior 
event. In the absence of events vetoed by the active anti- 
coincidence shield, this column would accurately reflect 
the true elapsed live-time, but the standard operating 
mode does not downlink vetoed events, so that in general 
adding the PRIOR column would not yield the proper 
live-time. However, for source count-rates significantly 
higher than the veto rate the error becomes negligible, 
and for the Crab we determine it is only ^0.3%, by com¬ 
paring the summed PRIOR column to the live-time re¬ 
ported by the instrument once per second. 

We used the live-time curve to adjust the net exposure 
time of each phase, and for each phase bin we extracted 
events for both 17" and 50" circular regions centered on 
the pulsar. The reason for picking two extraction re¬ 
gions is to test the stability of the fitting procedure. The 
nebular component will change as a function of the ex¬ 
traction region size, while the pulsar component should 
not, providing a cross check on how well the two compo¬ 
nents are distinguished. To avoid the statistical pit-falls 
of subtracting a large component, we decided to fit the 
two components together, presuming that each phase bin 
can be decomposed into the un-pulsed nebula, constant 
throughout all phases, and the pulsar. 

We first investigate what spectral models to use by fit¬ 
ting the three phase bins 10-12 corresponding to the 
pulse-off interval. We found that a power-law yields an 
inadequate fit (x^ed = degrees of freedom 

(dof)) and that a broken power-law provides a much bet¬ 
ter fit (x?ed = dof). We also tried out a 

logpar model and an exponential cutoff model, but 
both failed to fit the high energy tail of the spectrum 
and under-predicted the flux. Although we recognize 
that it is difficult to infer physical meaning from a broken 
power-law, we chose to continue with this representation 
because it gives a simple and intuitive understanding of 
the shape of the spectrum. 

We use the best fit broken power-law to investigate the 
shape of the pulsar component during the peaks in phase 
bins 2 and 8. Freezing the broken power-law model to 
what we derived for phase bins 10 - 12, we fit the pulsar 
component with a power-law and a broken power-law re¬ 
spectively. Once again the broken power-law clearly pro¬ 
vides a better fit; for phase bin 2 the power-law fit has 
Xred~^'^^ for 806 dof and broken power-law x?ed~^-^^ 
for 809 dof. On this basis we chose to model with two in¬ 
dependent broken power-laws. We fit a total of 26 spec¬ 
tra (13 phase bins for FPMA and FPMB) in XSPEC 
between 3-78 keV. 

We observe degeneracies between the models when we 
do not constrain the normalization of the pulsar compo¬ 
nent during phase bins 10-12. We resolve this by limit¬ 
ing the contribution of the pulsar during these phases to 
be a small fraction dur ing pulse-on phases. According to 


Weisskopf et al. (2011), the ratio of the 0.3 - 3.8 keV flux 
between the primary pulse peak and pulse-off is about a 


factor of 100, and we applied this restriction to the rela¬ 
tive pulsar normalization in the model for these phase 
bins. We coupled the break energy, BE, between all 
phase bins for the pulsar and nebula respectively. Table 

shows the best-fit parameters for the pulsar spectrum 
for both extraction regions in the top half for all phase 
bins except 10 - 12. The parameters for the two extrac¬ 
tion regions agree for most bins within Icr and all of them 
within 2cr. The bottom half lists the best fit parameters 
for the nebula in the two extraction regions and as antici¬ 
pated they are different, softening with increasing radius 
but maintaining the same break energy. Eigurej^ bottom 
panel shows the averaged power-law for each phase bin of 
the two extraction regions for Fi and F 2 respectively. We 
have omitted the parameters for phase bins 10-12 since 
the pulsar component could not be properly constrained 
during these periods. The average AF = r 2 — Fi across 
the phase bins, excluding bins 10 - 12, is 0.27 ± 0.09. 

3.3.2. Spatially Integrated Phase Resolved Speetrum 

We measured t he ph ase-averaged spectrum from a 200" 
region in Section |3.1| and found F = 2.0963 ± 0.0004. At 
smaller radii we find that the pulsar and the nebula can 
both be characterized independently by broken power- 
laws. It may seem strange that the superposition of ap¬ 
parent broken and un-broken power-laws should sum up 
to a simple power-law. To show how the Crab achieves 
this, we investigate here the phase-resolved total spec¬ 
trum (nebula-fipulsar) using the same phase bins as de¬ 
fined in the previous section. We use a broken power-law 
where required by the data, and power-law otherwise and 
the fit results are listed in Tabled 

Figure shows the ratio of the data to the phase- 
average power-law fit. The curves have been scaled 
for clarity to show how the Crab roughly decomposes 
into four components; 1st pulse, bridge emission, 2nd 
pulse, pulse-off (nebula). Although both pulse peaks 
have spectra that steepen (r 2 ^ 2.047 — 2.099), the 
spectra remain harder than the phase-aver aged value. 
The bridge emission is best approximated by a power- 
law with an index close to the phase-averaged value 
(F = 2.086T0.002), while the off-pulse nebula emission is 
significantly softer and mildly broken (Fi = 2.123T0.003, 
r 2 = 2.134(—0.004/ + 0.01)). The sum of nebula and 
pulsed emission thus approximately conspire to mimic 
the phase-averaged power-law and explain why as a 
whole we observe a power-law spectrum. 

Finally in Figure we show the spectral energy distri¬ 
bution for the (1) phase averaged pulsar + nebula, (2) 
nebula (phase bins 10 - 12), (3) pulsar + nebula (phase 
bins 1-9, 13), (4) and pulsar (phase bins 1-9, 13) 
with nebula subtracted. The fit parameters for (1) is 
F = 2.0963 ±0.0004, (2) F = 2.0865 ± 0.0004, (3) nebula 
alone as recorded in Table (4) and for the pulsar alone 
Fi = 1.72 ± 0.02, r 2 = 1.7T± 0.02 and BE = 10.5 ± 1.6. 

3.4. Spatially Resolved Speetroseopy of the Nebula 

Erom the phase-averaged spectra shown in Eigure 
it is evident the Crab spectral index changes quite dra¬ 
matically across the face of the remnant. This is in part 
due to the harder pulsar spectrum mixing with the softer 
nebula. When restricting only to phase bins 10-12 the 
change is less dramatic, going from ri/r 2 = 1.92/2.00 
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TABLE 2 

Phase resolved fits to extraction regions 17 " and 50" 

Tbabs(bkn+bkn) Tbabs(bkn+bkn) 

17": 15884 degrees of freedom) 50": 23752 degrees of freedom) 

Pulsar Pulsar 

bknpower bknpower 


Phase 

Pi 

r2 

BETkeV) 

Pi 

r2 

BE*' (keV) 

1 (0-0.07) 

1.64T0.04 

2.01T0.02 

11.7T0.6 

1.66T0.04 

2.01T0.02 

13.1T0.4 

2 (0.07-0.14) 

1.80T0.01 

2.02T0.02 


1.83T0.01 

2.05T0.02 


3 (0.14-0.21) 

1.52T0.04 

1.76T0.05 


1.55T0.04 

1.77T0.05 


4 (0.21-0.28) 

1.39T0.06 

1.62T0.07 


1.28T0.06 

1.64T0.07 


5 (0.28-0.35) 

1.30T0.05 

1.70T0.06 


1.29T0.05 

1.74T0.06 


6 (0.35-0.42) 

1.42T0.03 

1.75T0.04 


1.42T0.03 

1.76T0.04 


7 (0.42-0.49) 

1.57T0.02 

1.93T0.03 


1.61T0.02 

1.85T0.03 


8 (0.49-0.56) 

1.66T0.02 

1.92T0.03 


1.71T0.02 

1.91T0.03 


9 (0.56-0.63) 

1.66T0.08 

1.97T0.12 


1.83T0.08 

1.95T0.12 


13 (0.9-1.0) 

' °-0.08 

1.81T0.10 


1.96T0.08 

2.11T0.10 



Nebula Nebula 

bknpower bknpower 


Phase 

Pi 

r2 

BE (keV) 

Pi 

r2 

BE (keV) 

lO-lB (0.63-0.9) 

1.92T0.01 

2.00T0.01 

8.3T0.5 

1.99T0.01 

2.09T0.01 

8.3T0.2 


^Break energy (BE) component coupled for all phase bins. 
^Nebula component only. 



Energy (keV) 


Eig. 4.— Ratio plot of best fit model P = 2.1 of the 13 phase 
bins extracted from a region of 200" centered on the pulsar. 



Eig. 5.— Spectral energy distribution of the pulsar + nebula over 
all phase bins, nebula (phase bins 10 - 12), pulsar + nebula (phase 
bins 1-9, 13), and pulsar (phase bins 1-9, 13) alone with nebula 
subtracted. 


TABLE 3 

Phase resolved Pulsar+Nebula fits to 200" extraction 

REGION. 


Phase 

Pi 


pa 

BE^(keV) 

1 (0-0.07) 

2.090 

±0.003 

2.110 ±0.006 

11 ±3 

2 (0.07-0.14) 

2.041 

±0.003 

2 099“^°‘°®® 
^•'^^^-0.01 

12.3 ±1.2 

3 5-14-0.21) 

2.081 

±0.002 

- 

- 

4 (0.21-0.28) 

2.090 

±0.002 

- 

- 

5 (0.28-0.35) 

2.086 

±0.002 

- 

- 

6 5-35-0.42) 

2.064 

±0.002 

- 

- 

7 5-42-0.49) 

2.030 

±0.002 

2 Q47+0-006 

z.un/_o 01 

13.7 ±3 

8 5-49-0.56) 

2.034 

±0.003 

2.056 ±0.005 

10 ±1.5 

9 (0.56-0.63) 

2.120 

±0.002 

- 

- 

10 (0.63-0.7) 

2.120 

±0.004 

2.138l^;g'j'"‘ 

10 ±4 

11 (0.7-0.8) 

2.123 

±0.003 

Z.i9T_o 01 

10 ±4 

12 5-8-0.9) 

2.126 

±0.002 

- 

- 

13 5-9-1.0) 

2.124 

±0.002 

- 

- 

If r 2 and break energy (BE) 

is not given 

the best fit was a 


power-law. 


at 17"jl.99/2.09 at 50" (Table [|) to 2.12/2.14 at 200" 
(Table0. Softening of the spectra with increasing radius 
is p redicted in theory and has bee n observed in G21 .5- 


0.9 (jNynka et al.||2014|) and 3C 58 (iSlane et al. |2QQ4|) 
well in the Crab at lower energ:ies by Chandra JVVeissKi 


as 


eisskopf 
et al.|2QQQ ) and XMM-Newton ( Kirsch et al.|2006| . Due 


to mixing of spectra at different radial locations by the 
PSF in NuSTAR^ it is not straight forward to measure 
the true (unmixed) index as a function of radius, and a 
more careful analysis is required to confirm and quantify 
the effect. 

The theoretical predictions apply to the nebula, and we 
chose to analyze the spatial variation during the phase 
bins 10-12 when the pulsar is off. We slide a box of 
24.5" X 24.5" solid angle across the inner 100" to obtain 
a high S/N spectra, and the box is incremented along 
RA and Dec by a solid angle of 2.45", the size of a NuS- 
TAR projected sky pixel. At each step we calculate the 
responses for the box center and fit a broken power-law 
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Fig. 6. — Spectrum for off-pulse bins 10 - 12 extracted from map 
location (Ra,dec)=(0,0) of FPMA, which coincides with the pulsar 
location. The black and red data set is the same, but the models 
different. The red curve shows the best fit broken power-law, which 
had Fi = 1.91T0.01, F 2 = 2.03T0.01, E^reak = 8.7±0.9keV. The 
black curve is a power-law with index F = F 2 and the normalization 
scaled such that the high energy part is the same. 


spectrum. We scale the background from the simulated 
master background. We performed fitting using XSPEC 
and Cash statistics. 

In general a broken power-law adequately describes 
the data. However, in some cases the fitting found ei¬ 
ther Fi = r 2 , or a break energy less than 5keV. In 
these cases we reverted to a single power-law model. At 
the edges where the torus transitions into the extended 
nebula, the spectra become complex and can no longer 
be represented by a broken power-law either. Some of 
them show negative breaks and for the edge regions we 
therefore characterized the spectrum with a simple single 
power-law setting F = r 2 . 

In the following we switch to a coordinate system rela¬ 
tive to the pulsar location (ARa,ADec) = (0,0). To illus¬ 
trate the data quality and the level of deviation from 
a power-law in the center of the remnant. Figure 
shows the ratio of the spectrum extracted from FFMW 
for one grid point (24.5" x 24.5") with the center at 
(ARa,ADec)=(0,0) to the model. The red curve shows 
the best fit broken power-law, for which we find Ti = 
1.91 ±0.01, r 2 = 2.03 ±0.01, and E^reak = 8.7±0.9keV. 
The black curve shows a single power-law where we have 
fixed the index T = r 2 = 2.03 ±0.01 and normalization 
such that the spectra match at F> 10 keV. 

Simply by eye it is clear that these spectra strongly 
deviate from a power-law by softening above lOkeV. We 
emphasize that the effect of PSF mixing is to obscure the 
true underlying spectra. It cannot create the appearance 
of a broken power-law out of a superposition of spectra 
if they were solely power-laws. 

We present now the spatially resolved maps. Figure 
and Figure!^ show maps of the (1) the break energy, 

AF = r 2 — ii), (3) r<6/cey (the photon index at energies 
< 6keV), and (4) r 2 . In general the errors in individual 
boxes are (ri,r 2 ) ^ ±0.015 and the error in the break 
energy (if there is one) ±1 keV. The left panels show the 
map overlaid with NuSTAR intensity contours and right 
panel with Chandra contours to show details of the spa¬ 
tial structures. The pulsar is represented by a cross and 
located in all images at (ARa,ADec) = (0,0). Although 


we show the map for FFMA for ease of presentation, the 
maps from FFMB are identical within errors. 

Figure shows that inside the remnant the break en¬ 
ergy is roughly constant, with an average value of ^9 keV 
decreasing towards the NW and increasing at the SE 
edge. The AT map shows that the highest values of AT 
follows the curvature of the forward edge of the torus. 
It is interesting that the largest value is not found at 
the location of the highest intensity, but rather slightly 
more north on the edge of the Chandra intensity contour. 
Comparing to the break energy map, where AT is high 
the break energy is on average lower. 

Figure!^ shows the low-energy spectral index, r<6/cey- 
The 6 keV energy was chosen because it is always below 
the break energy. The greatest change in morphology 
in this map occurs along the forward edge of the torus 
where the spectrum softens rapidly. Along a line from 
the pulsar towards the NW corner, the radial spectrum 
softens faster above the break than below. This holds 
at the majority of azimuthal angles around the torus, 
but not in the jet direction. Neither the jet region or 
counter-jet region have measurable spectral breaks. The 
SE corner appears to have a negative AF, but this is 
likely the effect of the FSF scattering the harder torus 
high energy spectrum into the softer nebula. 

There is no easy way to spatially disentangle spectral 
components that have been mixed by the FSF (In the 
Appendix we provide details on the NuSTAR FSF and 
discuss the technical details of it). To interpret the spec¬ 
tral map, we therefore compose a model of the nebula 
and forward fold it through the NuSTAR response and 
compare it to the data. We create a 2-D sp ectral model 
based on the analysis of |Mori et al.| (|2004| of the Crab 
with Chandra. This data set shows spectral variations on 
arcsecond scales, but we follow the authors comment that 
the remnant can be represented by six components: (1) a 
halo (shown in dark red in Figure 1^ with T = 3.0, (2) a 
cap (yellow) T = 2.5, (3) a skirt (li^t blue) T = 2.1, (4) 
a torus (dark blue) T = 1.9, (5) a jet (medium blue) 
T = 2.0, and (6) a center (black) T = 1.6. We as¬ 
sume that this spectral model holds up until our mea¬ 
sured break. Based on the NuSTAR analysis, we fur¬ 
ther assume that the torus has an average spectral break 
at 9keV with a AF = 0.25, and that the spectrum of 
the center has a break at the same energy, but with a 
AF = 0.1. We arrived at this smaller AF = 0.1 through 
a series of iterations since this component is maximally 
obscured and difficult to determine from the NuSTAR 
data. The center and the torus are the only two compo¬ 
nents that we allow to have a break. All other regions are 
presumed to have power-law spectra. We used a Chandra 
image of the Crab, shown in Figureto set the normal¬ 
ization of the map in the 2 - lOkeV band. The spectra 
are completely defined by this set of parameters. 

We propagated this model through the NuSTAR FSF 
and responses and analyzed the output in exactly the 
same manner as the real data itself. Figure shows 
the resulting maps. Macroscopically they reproduce the 
NuSTAR observations very well. There are some discrep¬ 
ancies, but these can be explained by the simple model 
we use whic h does not i nclude the low energy variations 
observed by Mori et al. (2004) or the break energy varia¬ 
tions seen by NuSTAR. However, these simulations still 
address some important questions. It is clear that the 
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Fig. 7. — Plots showing the break energy and AF map of F 2 -F 1 for FPMB during pulse off. Left panel contours are the NuSTAR intensity 
levels and the cross marks the pulsar location. Right panel contours are from Chandra. 


torus spectrum has to steepen with energy, and that one 
has to use an approximately constant spectral index in 
both Fi and r 2 all the way out to the edge, where there 
is a rapid transition to a softer, unbroken spectra. We 
tested this hypothesis by creating a model with a lin¬ 
early changing spectral index as a function of radius as 
has been o bserved for G21.5-0 .9 (Safi-Harb et al. 2001) 
and 3C 58 (Slane et al. 2004). 'I'he resulting maps do 
not match the observations and create spectral slopes 
that are far too soft in the interior. The spectrum of the 
core m ust be hard i n ord er to reproduce the NuSTAR 
results. |Mori et al.| (|2004|) did not find a hard center of 
the remnant below 10 keV, but their observations suf¬ 
fered pile-up in the center, and they couldn’t measure 
the spectrum in the inner 10". The origin of the hard 
component seen by NuS TAR could be the neutro n star 
itself, which according to 


spectral index of F ' 


Weisskopf et al. (2011) has a 
1.9 ±0.4. Because the simulations 


are driven by the Chandra flux map, the location of the 
pulsar has no counts because it was excised, which may 
explain why instead of a point source our simulations re¬ 
quired an extended, hard core, distributing the missing 
flux over a large region. Finally neither the cap or the 
jet appear to steepen with energy 

It is pertinent to question whether dust scattering 
could have anything to do with the slope below 10 keV. 
Unless the Nh column inside the remnant is larger than 
the galactic column by a factor of 10, the interstellar dust 
extinction has little effect above 3keV (Seward, Goren- 
stein fc Smith|[2QQ6 ). Indeed the fact that the break en¬ 
ergy and Al' trace out features of the remnant strongly 
indicates that it must be intrinsic to the source. 

3.5. Spatial Extent of the Nebula 

The size of the Crab remnant shrinks as a function 
of energy due to the radiative lifetimes of outward- 
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Fig. 8. — Plots power-law map of r<6fcel/ and r2 for FPMB during pulse off. Left panel contours are the NuSTAR intensity levels and 
the cross marks the pulsar location. Right panel contours are from Chandra. 


propagating electrons being shorter for high energy than 
low energy particles. This effect is often referred to as 
‘synchrotron burn-off’. To investigate the radial extent of 
the Crab as a function of energy we deconvolved the NuS¬ 
TAR maps using a maximum likelihood method. The 
deconvolution procedure is sensitive to artifacts, such as 
detector gaps and the variation in signal to noise ratio 
with position in the map. The stronger the source rela¬ 
tive to background the better the deconvolution results. 
The PSF is relatively constant near the optical axis (the 
difference in the HPD between off-axis angles of V and 
2' is less than 1"), but becomes azimuthally distorted 
at large off-axis angles. This, however, does not become 
noticeable until about 3' off-axis, where the difference be¬ 
tween the major and minor axis of the PSF is ^2%. To 
minimize these effects we selected a subset of the obser¬ 
vations in Table [^marked with “b” at off-axis angles less 
than 2' and well away from the detector gaps. At these 


off-axis angles we can design an average PSF weighted 
by time and combine the images to yield a more robust 
result than individually deconvolving short segments. 

To remove the contamination of the pulsar we only use 
photons falling in phase bins 10-12 and deconvolve the 
Crab in the following energy bands: 3-5,5-6,6-8,8- 
12, 12 - 20, 20 - 35, and 35 - 78keV. Prior to combining 
the images, we vignetting correct them with the effective 
area taken at the area-weighted average energy. 

Deconvolution with a maximum likelihood method is 
iterative and if not performed with care can introduce 
artifacts due to over deconvolution. During the deconvo¬ 
lution process we checked the relative size between the 
selected iteration steps (20, 30, 40 and 50) within each 
energy band, and saw that the size difference remained 
largely unchanged after 30 deconvolutions. For the high¬ 
est energy band, however, we started to observe artifacts 
after 50 deconvolutions, and we therefore estimate that 
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Fig. 9. — Top: Input power-law index map for the simulations. 
For more details see description in text. Bottom: Chandra inten¬ 
sity map used to set the 2 - lOkeV flux level for the simulations. 


40 deconvolutions is a safe number to use for all energy 
bands. 

Because of the nature of the deconvolution process, it 
is difficult to assign an error to the resulting size. Ideally 
a model of the source can be forward folded and com¬ 
pared to the actual image, but in the case of the Crab 
where such a model is not well known, and the morphol¬ 
ogy of the source is complex, this approach is not fea¬ 
sible. Forward folding the deconvolved image with the 
PSF and comparing to the raw image is another option, 
but it only serves to identify gross errors. It is therefore 
not possible to assign an error to the absolute size in 
the deconvolved images. Fortunately, we are interested 
in a rate of change as a function of energy, and we can 
investigated the error of the deconvolution as a function 
of energy by deconvolving several strong NuSTAR point 
sources. These should have a constant radial extent as 


TABLE 4 
Nebula HWHM 



FPMA 


FPMB 


Axis 

7 

k 

7 

k 

NE Torus 

0.094A0.018 

49(1) 

0.079T0.017 

50(1) 

SW Torus 

0.060A0.020 

38(1) 

0.085T0.02 

42(1) 

SE Jet 

0.083T0.062 

13(1) 

0.014A0.056 

12(1) 

NW Jet 

0.245A0.029 

46(2) 

0.192 ±0.027 

42(2) 


a function of energy, and any discrepancy is assumed to 
be the error introduced by deconvolution. We used the 
same energy bands as for the Crab images, and found an 
evolution in the deconvoled point source sizes as a func¬ 
tion of energy. The discrepancy between the highest and 
lowest energy bands is ^1.5", and we have conservatively 
assumed a 2 " la error in the relative size of the remnant 
between the energy bands. 

Figure \TT\ shows the resulting images for FPMA. The 
major components observed in Chandra are also seen by 
NuSTAR, and morphological changes as a function of 
energy are clearly visible. Figure plots the contours 
of the HWHM as defined by the point where the inten¬ 
sity falls to half the value measured at the pulsar loca¬ 
tion, and shows the magnitude of the shrinkage is posi¬ 
tion angle (PA) dependent. Figure 13 shows the profiles 
of two perpendicular slices with the profiles extracted 
by bi-linear interpolation. The intensity map is off-set 
from the pulsar due to the 27° torus inclination and the 
beaming of the leading edge of the torus, but we nor¬ 
malized the curves at the pulsar location, which places 
the peak intensity towards the NW for low energies. The 
offset decreases with increasing energy and is gone by 20 
keV. The NE and SW side of the torus plane therefor 
have somewhat different morphologies close to the pul¬ 
sar, and it can also be seen that the jet and counter-jet 
directions evolve quite differently; the counter-jet and 
cap area HWHM falls off rapidly with energy, while the 
forward jet is much slower. 

We fitted the HWHM of the profiles as a function of 
energy using a power-law kE~^. Figure 14 shows the 
fits for FPMA for the two sides of both the jet axis and 
torus plane, and Table summarizes the fit values for 
both FPMA and FPMB. The flux of the jet is much 
smaller than that of the torus, and what we measure 
along the SE is most likely not the jet, but rather the edge 
of the torus. The rate of change in NE, SE and SW are 
of a similar slope, but different intensity, indicating the 
shrinkage occurs in the plane of the torus. The projection 
effect of the tilted plane would in this way explain the 
smaller magnitude observed along the SE edge. The NW 
edge of the torus has a line of sight that also includes the 
NW cap and counter-jet, obscuring the true shrinkage of 
the torus in this direction. We know this area has the 
softest spectrum and the index rate of shrinkage, 7 , is 
twice as large as that of the torus. 


4. DISCUSSION AND CONCLUSION 


Gratton (1972) and Wilson (1972) developed the first 
analytic niodels attempting to explain the observational 
properties of the Crab. They constructed a model of the 
nebula assuming that electrons and positrons produced 
by the pulsar are transported away by diffusion, with 
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Fig. 10 . — Simulated model. Left panel: Break energy and AF map of F2-F1. Right panel: power-law map of F<6fcel4 and F 2 . Contours 
are from Chandra 


energy lost through synchrotron radiation. The size and 
spectral shape of the nebula in the radio and optical are 
well explained by the model, but the authors did not 
consider the details of the production of the electrons 
and positrons by the pulsar, or their relationship to the 
p ulsar properties, such as i ts magnetic field. 

Rees & Gunn (1974) and [Kennel fc Coroniti ( 1984a|b ) 


constructed a steady state spherically symmetric mag- 
netohydrodynamic (MHD) model of the Crab with a 
toroidal magnetic field, which links the pulsar to the neb¬ 
ula. A highly relativistic pulsar wind is terminated by 
a strong MHD shock, which decelerates the flow to non- 
relativistic speeds. Downstream the approximately adi¬ 
abatic flow continues to decelerate, carrying its frozen- 
in magnetic field out to the supernova ejecta that con¬ 
fines the synchrotron nebula. The model successfully de¬ 
scribes the integrated spectrum from optical to X-rays 
where the diffusion model fails. However, it does not re¬ 


produce the radio spectrum, nor reproduce the spatial 
variations of the optical spectrum across the remnant. 
The model predicts a constant spectrum that steepens 
only at the edge of the remnant, but this has not been 
met by observations in the optical, which instead show 
a monotonically increasing spectral i ndex (Veron-Cetty 


& Woltjer 1993 Temim et al. 2006). Siniilar discrep 
ancies have been observed in A-ray s for the other two 


young PWN, 3C 58 and G21.5-0.9 (jNy nka et al. 2014 
Slaneet^2l22M _ Safi-Harb et al.||200l| ). 4’his motivatec. 


'I'ang Chevalier 2012|) to re-inve stiga te the d i ffusion 


driven models from Gratton (1972) and Wilson (1972). 
They found that a combination of diffusion and advection 
can reproduce the observed spectral index variations ob¬ 
served in radio and optical for the Grab, and in X-rays for 
the two other sources. However, diffusion under-predicts 
the size of the Crab in X-rays. A combination of ad¬ 
vection and diffusion increases the predicted half-light 
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Fig. 11.— Maximum likelihood deconvolved images of FPMA in six different energy bands shown in a square root stretch during pulse 
off. The shrinking of the NW counter-jet is easily observable and the morphological changes of the NE torus can also be seen. 



Fig. 12.— Image of the HWHM contours for energy bands: 3 - 
5, 5 - 6, 6 - 8, 8 - 12, 12 - 20, 20 - 35, and 35 - 78keV. 


radius of the Crab at 10 keV from 20 " in the pure diffu¬ 
sion case to 25" for the best fit combination of advection 
an d diffusion. This is still t oo small, and the conclusion 
of Tang & Chevalier ( 2012 ) is that due to the compact 
size of the Crab in X-rays, the advection time scale dom¬ 
inates over the diffusion time scale for X-ray producing 


particles. Detailed Chandra maps of the Crab nebula by 


Mori et ah (2004) have shown that the spectral index, 
within some small-scale variations, is constant for most 
of the nebula and steepens abruptly at the edge. This 
supports the idea of an advection, rather than diffusion, 
dominated X-ray nebula. 

The first measurements of the extent of the X-ray Crab 
as a function of energy from 2-12 keV were made in 


1964). 


Ku et al. ( 

19761) later combined the results of 

several experiments (iFalmieri et al.||1975l 

Ricker et al. 

1975 

Fukada et al. 1976 Ku et al.| 1976 

), and found 


effect they attributed to a combination of diffusive and 
advective transport of electrons in the nebula. At this 
time the advective model of KC84 had not yet been con¬ 
structed, but KC84 concluded in their paper that these 
results agree fairly well with an advective electron trans¬ 
port. These rates were obtained from multiple experi¬ 
ments and observatories along a narrow range of posi¬ 
tion angles through the Crab, and they did not probe 
any azimuthally dependent information. 

With NuSTAR we have measured the energy depen¬ 
dent size at all position angles at energies from 3 - 
78keV. The spatial dependence of the HWHM along the 
edge of the torus plane is well fit by a power-law with an 
average of 7 = 0.08 ± 0.03. This is in good ag r eement 
with the predictions made by Kennel & Coroniti (1984b) 
of 7 ^ 1/9. Measurement of the jet rate of change could 
not be done since the torus flux dominates, but judging 
from the tail of the profile in Figure bottom panel, it 
does appear as if the jet shrinks more swiftly outside the 
torus as a function of energy. The size of the counter-jet 
as a function of energy clearly follows a different rate. 
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Fig. 14. — Half Width Half Max, HWHM, in arcseconds. The 
HWHM is measured from the pulsar location along the torus and 
jet axis for either side. From top to bottom: NE torus plane. SW 
torus plane. NW jet axis. SE jet axis. 



Jet axis (") 


Fig. 13.— Energy dependent profiles of module A interpolated 
along the torus plane and jet axis of the deconvolved images shown 
in Figurein the direction of the arrows. Top: Profile along torus 
plane. The right shoulder corresponds to the SW. Bottom: Profile 
along jet axis. The right shoulder corresponds to the NW counter¬ 
jet and shows a very clear shrinking as a function of energy, while 
the left shoulder corresponds to the SE jet and shows almost none 
at all. 


7 = 0.22±0.03, suggesting the energy transport is not 
the same as in the torus. In other band passes this re¬ 
gion is found to be different as well, and is the only part 
of the remnant where the synchotron nebula extends be¬ 
yond the edge of the visible filaments. It is purported 
that here the shock velocity is larger than for the rest 
of the remnant, and the post-shock cooling time longer 
than the age of the Crab ( Sankrit fc Hester]|1997 ). 

Similar values of 7 ^ 0 . 2 , h ave been found for tne neb¬ 
ula remnants in MSH 1 5-52 (|An et al.||2014|) and G21.5- 
0.9 (Nynka et al. 2014), two FWIN observed with NuS- 
TAn and covering the same energy band. The authors 
conclude that both remnants appear to be diffusion dom¬ 
inated. Both of these remnants are more than twice as 


distant (^5 kpc) and have distinct morphological differ¬ 
ences from the Crab. MSH 15-52 is not a spherically 
symmetric remnant, but highly irregular and dominated 
by a powerful jet. G21.5-0.9 is spherically symmetric, but 
the inner portions of the remnant are not resolved and 
there is no apparent torus or jet evident. This would sug¬ 
gest that for the Crab we are probing the inner workings 
of the PWN, namely the torus where advection processes 
dominate, while in the counter-jet region diffusion pro¬ 
cesses dominate similar to what is found in MSH 15-52 
and G21.5-0.9. 

Using a map of the average spectral index distribu¬ 
tion from Chandra as the basis, we have shown through 
simulations that the torus must have a near constant 
spectrum as a function of radius. If we presume that 
the wind expands radially along the torus plane, then 
the transition from a flat profile to a steep occurs at a 
radius of approximately at 50". This is consistent with 
the KG84 model, which predicts a flat profile in an ad¬ 
vection dominated MHD flow that transitions at ^ 50" if 
the termination shock is of size 10". In this geometry, the 
energy dependent rate of shrinkage we have found is also 
consistent with KG84. We confirm that the torus is well- 
described by a portion of a spherical outflow, with rel¬ 
atively little diffusion or turbulent mixing and conclude 
that for the X-ray band at least, advection dominates the 
processes. 

What was not anticipated is that in the NuSTAR band 
there is a very clear spectral steepening inside the rem¬ 
nant. Our analysis has revealed the spatially resolved 
nebula is not well fit by a power-law across the band 3 - 
78 keV, but rather needs a model with a spectral steep¬ 
ening. We chose a broken power-law as an adequate rep¬ 
resentation and presented maps of the spectral variation 
that show a steepening of AT ^ 0.25, which appears to 
be limited to the torus of the nebula with an approx¬ 
imate break energy at ^ 9keV. We confirm this with 
simulations in order to address any ambiguities resulting 
from spatial smearing by the PSF. We also conclude that 
the jet and counter-jet regions do not appear to have a 
similar spectral steepening. 

The integrated spectrum of the nebula has been known 
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to turn over above 100 keV, transitioning from a F ^ 2.1 
to ^ 2.14. The gradual break must come about from the 
shrinking of the nebula, and the steepening of the torus 
spectrum. There is no immediate physical interpretation 
of the torus spectrum. It is possible the steepening could 
result from the projections of the emission from electron 
populations of different synchrotron ages, but a similar 
steepening has also been by NuS TAR in the spatially 
integrated spectrum of G21.5-0.9 (|Nynka et al.||2014|), 
and since the two remnants are morphologically diherent, 
it questions the geometrical argument and might suggest 
the steepening is connected to the injection spectrum 
itself. In either case the steepening of the torus spectrum 
shouldn’t come as a surprise. Hard X-ray instruments 
have measure d the photon index above 100 keV to be 
2.140 ± 0.001 dPravdo, Angelini & Harding||1997|) while 
Chandra me asures the average photon index of the torus 
to be ^1.9 (Mori et al. 2004). This is a softening of 
AF ^ 0.25, which matches nicely with the number found 
in our maps. Using the rate of burn-off for the NE side 
from Table [4j the HWHM radius of the Crab is just ^30" 
at 100 keV, restricting the source of the steepening to 
come from the innermost regions of the remnant. It is 
therefore not only likely but necessary for the torus to 
steepen in the NuSTAR band in order to bridge the gap 
between soft X-ray and y-ray observation. 

We performed phase resolved spectroscopy of the Crab 
on several different length scales from the inner 17" out 
to 200". The pulsed spectrum is best represented by a 
steepened spectrum with a break energy of ^ lOke V. As 
found previously b y 
ingale et al. ( 2001 ), 
spectral evo 


Weisskopf et al. (2011) and Will 
he index below lOkeV, Fi, shows 
ution as a function of phase; the secondary 
pulse is harder than the primary and the hardest index 
occurs during the bridge emission between the first and 
secondary pulse. We find that F 2 traces Fi with an av¬ 
erage AF = 0.27 ± 0 .09. There are indications of this 
steepening in RXTE (Pravdo, Angelini & Harding 


observations , but it was not properly quantified. [Kniper 


et al. ( 2001 ) present phase resolved pulsed spectra by 


combining BeppoSAX^ CGRO^ and GRIS data and found 
that from 0.1 keV up to 10 CeV it could be fit using three 
components: ( 1 ) a power-law, ( 2 ) a modified power-law 
(logpar) for the first pulse, and (3) a modified power- 
law for the bridge emission. They found acceptable fits 
by combining these three models, and we attempted to 
fit with the same combination. While we were able to 
find statistically acceptable fits, the models proved de¬ 
generate in our narrower band-pass without the 7 -ray 
spectrum to constrain them. 

Einally NuSTAR participated in a ToO observation of 
the Crab during the flaring on the 9th of March 2013. 
No flux change was detected between 3-78 keV aside 
from what is expected from calibration, which is ±5% in 
flux and 0.01 in spectral index change, and thus places 
an upper limit on the hard X-ray variability due to 7 -ray 
flares. 


This work was supported under NASA Contract No. 
NNC08ED60C, and made use of data from the NuSTAR 
mission, a project led by the California Institute of Tech¬ 


nology, managed by the Jet Propulsion Laboratory, and 



Fig. 15.— NuSTAR encircled energy curve at 10 keV. 

istration. We thank the NuSTAR Operations, Software 
and Calibration teams for support with the execution 
and analysis of these observations. This research has 
made use of the NuSTAR Data Analysis Software (NuS- 
TARDAS) jointly developed by the ASI Science Data 
Center (ASDC, Italy) and the California Institute of 
Technology (USA). 


5. APPENDIX 

The NuSTAR point spread function has a sharp core 
(EWHM = 18") but extended large wings. It has an 
energy dependency, which causes the half power diameter 
to shrink as a function of energy by a few arcseconds 
between 3 and 10 keV. Above 10 keV the PSE remains 
constant. 

As shown in Eignre p!5]of an encircled power curve of 
a 10 keV PSE, 30% of the photons of a point source are 
located within a radius of 20", 50% within 35" and 80% 
within 70". This means that for an extended source like 
the Crab, which has most of its flux contained within an 
ellipse of 100"x 150", the extended wings will roughly re¬ 
distribute 60% of the flux from a central volume element 
over the rest of remnant. This causes the true spectral 
distribution to be mixed, and even with detailed knowl¬ 
edge of the PSE and effective area, it is unfortunately not 
possible to disentangle these mixed spectral distributions 
by any method of deconvolution. 

One approach to tackle this mixing, is to subdivide the 
extended source into smaller regions and calculate the 
cross-correlation functions between t he regions of both 
the AREs and PSEs as was done in |Wik et al.| (|2014|). 
The resulting spectra and responses must then be fit- 
ted simultaneously. This works well for large sources 
with spectra that are slowly changing, but for for smaller 
sources with fast changing spectra where the number of 
regions could number in the hundreds, this starts get¬ 
ting cumbersome. Another way is to forward fold a two 
dimensional spectral model through the optics response, 
carefully taking into account the PSE, and compare to 
actual data. 

Based on the fine structure we observe see in the spec¬ 
tral data maps, we have in chosen the latter method. 
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